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NOMENCLATURE 


Injection Hole Diameter 
Enthalpy 


Heat Transfer Coefficient with Film Cooling 
(spanwise averaged), Q"/(T,, - T,) anc WT. TA 


Heat Transfer Coefficient with Film Cooling 


(spanwise averaged), GRACIE e - T) 


Heat Transfer Coefficient without Film Cooling 
= Z -ٌ 
A and q / رت‎ L 
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(spanwise averaged), GORE 
Injection Momentum Flux Ratio, له ده‎ 
Thermal Conductivity 

Mole Fractions of Constitutant Gases 
Blowing Ratio, M o 

Static Pressure 

Molecular Prandtl Number 

Wall Heat Flux 

Reynolds Number, ۵ ما لام‎ 

Absolute Temperature 

Mean Velocity 

Velocity Component in X Direction 
Velocity Component in Y Direction 


Molecular Weights of Constitutant Gases 


Distance from Downstream Edge of Injection Holes 


Nondimensional Distance 
Thermal Diffusivity 
Density 


Adiabatic Film Cooling Effectiveness (spanwise 


averaged), eqn. 15 and 22 

Absolute Viscosity 

Kinematic Viscosity 

Nondimensional Coolant Temperature, eqn. 16 


Injection Angle Measured from Mainstream 


Flow Direction 
Nondimensional Distance Parameter, eqn. 25 


Second Coefficient of Viscosity 


SUBSCRIPTS 


Adiabatic Wall 
Coolant 

Some Reference Value 
Total — 
Recovery Condition 


Freestream 


SUBERSERERTS 
- Averaged or Mean Value 
- Variable Property 


- Fluctuation from the Averaged Value 
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x/d 


aw 


I. INTRODUCTION 


Over the past thirty years, gas turbines have become 
smaller and lighter with improved fuel economy. These 
desirable characteristics have been achieved by raising the 
turbine entry temperatures. To prolong useful blade life in 
this high temperature enviornment, cooling systems, thermal 
barriers, and/or different materials capable of withstanding 
higher temperatures are utilizied. 

This thesis concentrates on film cooling. Eilm cooling 
protects components from high gas temperatures by blanketing 
the component with a film of cooler air. This film insu- 
lates the component and acts as a heat sink. The problem 
encountered with film cooling is the prediction of the 


component's wall temperature. Two different convective heat 


transfer coefficients are currently used for this 
prediction: 

a چاو‎ 2 | (1) 
and 

0۳ = wre - T (20) 


Here, different temperature differences are used to define h 
and hg. 

Ihe objective of this thesis is to find a relation 
between the two heat transfer coefficients when the effects 


of variable properties are important. 
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11. BACKGROUND 


e,  1DRDENE BEADE COOLING FUNDMENTALS 

For a given size gas turbine, raising the turbine entry 
۳۱۳ results in a higher specific power. This 
enables the designer to operate the ‘engine cycle’ at a 
higher pressure ratio, resulting in improved thermal effi- 
Ciency of the thermodynamic cycle and improved overall effi- 
CIENCY: The net result is a lower specific fuel 
consumption. 

Raising the turbine entry temperature results in higher 
component temperatures. These components become more prone 
to failure at higher temperatures through greater oxidation 
rates of their surfaces, excessive creep, and/or by thermal 
fatigue. Figure 1 (Hennecke, 1382) shows a comparsion of 
the rising turbine entry temperature to the rising blade 
material temperature limitations. Notice that uncooled 
metal components are limited to a turbine entry temperature 
of approximately 1200 K. Current fuel economy goals require 
a turbine ۶۰۶۶ temperature of 1600 K and greater. 
Numereous cooling systems have been devised to maintain 
internal components at or below their maximum allowable 
temperature. Vanes, blades, disks, liners, and casings may 
all require some form of cooling protection. 

Figure 2 (Daniels and Schultz, 1982) illustrates two 
possible cooling medium schemes. Liquid systems use 
external sprays and/or internal convection methods. Air 
schemes use external films and/or internal convection 
methods. Liquid systems have several thermodynamic advan- 
tages over air systems. Dispite these advantages, however, 
air systems have been universally accepted due to their 
relative mechanical simplicity to implement. Three prin- 


Ċipletmmethodsumor air cooling currently exist. They are 


JJ: 


internal convection; film and -transpiration /efI us Tomni 
Figure 3 (Hennecke, 1982) illustrates these characteristic 
types of cooling for turbine blades. 

Internal convection is cooling by convection between the 
internal sections of the blade and the air. This type of 
cooling may be subdivided into several methods consisting 
of: 


(315) کت‎ through radial holes or channels of a single 
Or multipass' nature. 


(2) Enhanced heat transfer due to an enlargement of the 
internal surface area, promotion of .turbulance, or 
impingement ; 


(3) An insert that, distributes the air according to the 
localized cooling requirements. 


Film cooling is localized ejection of air through 
discrete slots or holes on the blade. This results ina 
film of air around the blade which is cooler than combustion 
gas temperature, acting as a heat sink. The film also acts. 
like a thermal barrier between the hot gases and the blade's 
surface. This tvpe of cooling mav be subdivided into 
several methods consisting of: 


(1) Single films; where specific parts of the blade are 
covered by a film; 


(2) Full coverage; where many films are placed close 
enough together so they cover part or all of the 
lade. 

Transpiration or effusion cooling is the bleeding of air 
through the blade at every point to form a protective film. 
The blade is made of a porous material to accomplish this 
method of cooling. 

Typically, cooling air is extracted from various stages 
in the compressor and supplied to the components that 
require cooling. Figure 4 (Hennecke, 1982) illustrates the 
coolant flow path for a single stage of a typical gas 
turbine. Typical coolant flow rates as a percentage of the 
gas turbine's total flow rate are included. Notice that the 
guide vanes and turbine blades have the highest coolant flow 


rates. These two components are of primary concern for they 
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have the largest amount of surface area exposed to the main- 
stream flow and consequently the highest thermal loading. 
Current applications use a combination of the internal 
convection and film cooling methods. ۳ 1 ۳ SS (Daniels and 
Schultz 


010 17021٢ Cooling air comes in from the top and bottom 


eo.) mipliuctiatesma  cCypmed. cuilde vane cooling 
of the vane. Notice the internal convection methods used: 
صا‎ is distributed to the internal surface; 

(2) Impingement on the surface after distribution; 

(3) Ribs and pedestals used to enlarge the surface area. 
After the air has cooled the vane by internal convection, it 
is ejected into the mainstream and-used as a film on the 
exterlur surface. Initial stages utilize the internal 
Samvection and film cooling principles as illustrated in 
Faroe 5. Latter stages may employ internal convection 
cooling solely due to the lower gas temperatures encountered 
in these stages. 

The goal in using any cooling system is to maintain the 
internal components temperatures at or below their maximum 
allowable temperature. Achievement of this goal will result 
in high reliability and long service life. In addition, it 
is desirable to obtain (as closely as possible) a uniform 
blade temperature distribution to minimize thermal stresses. 

In the selection of a cooling system to be employed, two 
key limitations must be kept in mind. These limitations 
are: 

(1) The minimal amount. of cooling air should. be used to 
minimize it S detrimential effect on engine perform- 
ance; 

(2) The effectiveness of a particular cooling, method to 
keep the component sufriciently coole for the 
desired turbine entry temperature. 

Ihe detrimental effect on engine performance is two 
fold. First, extracting air from the compressor reduces the 
available output power. Some work can be recovered by rein- 


5-62 he cooling air back into -the mainstream flow 


ns 


20۶667601666 ling. Second, ejection of air through the blade 
has a detrimential effect on the aerodynamic efficiency of 
the turbine. Diminished aerodynamic efficiency results from 
pressure losses and flow separation in the wakes trailing 
the blades. 

Figure 6 (Hennecke, 1982) shows the effectiveness of 
different methods used to cool components as a function of 
the turbine entry temperature. "Cooling effectiveness' is 
defined as the ratio of actual to theoretical blade tempera- 
ture reduction. This 'cooling effectiveness' is not only a 
function of the geometric configuration involved, Sout aco 
the coolant flow rate. As the turbine entry temperature 
increases, so must the required cooling effectiveness. 

Figure 7 shows the cooling flow required by various 
cooling methods as a function of the turbine entry tempera- 
ture. Notice that convection cooling requires the most 
cooling flow rate while effusion requires the least at any 
given turbine entry temperature. Although this figure is 
for specific leading edge and coolant temperatures, the 
trend is the same for other temperature combinations. ll 

Figures 6 and 7 indicate that to achieve the highest 
possible turbine entry temperature with minimum detrimential 
effects, effusion cooling should be used. Effusion cooling, 
however, is not currently used due to structural and manu- 
facturing difficulties, and because the porous blades often 
become clogged. This is one reason why film cooling methods 
are currently employed. 

Logically, one could assume that the more coolant pumped 
into the film and the fuller the coverage, the higher the 
turbine entry temperature. However, there are limitations: 

(1) The detrimential aerodynamic effects of film cooling 
may offset the gains of a higher turbine entry 
temperature;j 


(2) "Lift-Off" may occur where the film ceases to protect 
the surface. 
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Thus, thermal design of vanes and blades with film cooling 


must be optimized. 


Be  FIEM GOOLING 

Film cooling is the ejection of coolant into the 
boundary layer on the blade. Figure 8 shows a typical film 
cooling geometry with temperature profiles for a row of 
holes. Notice on the top view, at some distance downstrean, 
the flows from adjacent jets converge to form a continuous 
film over the blade's surface. In this region, lateral 
variations become negligible, allowing a two dimensional 
analysis to be conducted with reasonable accuracy. 


The typical temperature profiles on Figure 8 illustrate 


the effect of film cooling on heat transfer. The left 
profile is without film cooling. Heat transfer is from the 
freestream to the blade's surface. The driving force is a 


function ‘of the freestream recovery temperature and wall 
temperature. The right profile is with film cooling. Heat 
transfer is from the freestream to the coolant and from the 
coolant to the blade's surface. The "driving force iS a 
function of the freestream recovery temperature, coolant 
recovery temperature, and wall temperature. As the flow 
moves downstream from the injection site, the coolant mixes 
with the existing boundary layer. Eventually a temperature 
profile similar to the left will prevail. When this 
happens, the heat sink and insulating effects of the film 
become minimal. 

Figure 9 illustrates a typical heat transfer coefficient 


and adiabatic wall temperature distribution around a turbine 


blade. Notice that the highest temperature and heat 
transfer coefficient occurs at the nose of the blade, near 
the stagnation point, and the location of the highest 
thermal loading. Because of this, injection sites are 


located in the nose and utilize the lowest coolant tempera- 


ture. In the region around the stagnation point, the heat 
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transfer coefficients in Figure 9 are representative of a 
laminar boundary layer. At approximately one-fourth the 
distance along the suction surface, transition to a EUGENE 
lent boundary layer begins, resulting in an increase in the 
adiabatic wall temperature and heat transfer coefficient 
approximately three-fourths the distance along the suction 
surface. The heat transfer coefficients along the pressure 
side in Figure 9 are representative of a turbulent boundary 
layer developing with a constant streamwise pressure 
gradient (after an initial drop from the stagnation point). 
Near the mid-chord position, the heat transfer coefficients 
begin to rise as a strong favorable pressure gradient 


develops. 


C. BOUNDARY LAYER EQUATIONS 


To predict the heat flux from the mainstream to the 


blade, the equations governing the physical situation must 
be solved. Figure 10 illustrates the coordinate system 
used. Taking a differential element and applying the 


conservation principles of mass, momentum, and energy; the 
governing equations are derived. Using the governing equa- 
tions, a solution to the heat transfer problem is considers 
ably difficult since flow and temperature distributions are 
very complex. These equations are significiently simplified 
by making the following assumptions: 


Cl)» Ju AV 


AP ㆍ 
(2) ZL 
dy 
aT oT 
(3) LAS 
dy 9X 
Û 
(4) Ã + == 0 
S 


(5) No internal heat sources exist 


(6) No chemical reactions occur within the element 
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(7) Body forces are negigible 
(8) No work is performed by external fields 
Applying these assumptions reduces the governing equa- 
tions to the boundary layer equations. For turbulent 
boundary layers, the following form is’ for time-averaged 


velocities and variable properties: 








CONTINUTTY 
9 _ 9 
d E ۵ (5) 
X dy 
MOMENTUM 
_ 66 sq SDA mo Se 
1 نا‎ ٤ لا‎ 이 — —(u-c- pu iv”) (4) 
ax dv dx dv ay 
ENERGY 
3H 3H _dp 3q 3W 
imi Nu + T— (5) 
ö X d dy ax dy dy 
where 
du 
T = y—- ou'v' (6) 
dy 
and 
3T 
q pe ov'H' (7) 
d y 
The boundary layer equations, as presented, are highly 
nonlinear. The momentum and energy equations are coupled 


through the temperature dependency of density and viscosity. 


Thermal conductivity is also temperature dependent. The 
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dependence of these transport properties on temperature 
results in different velocity/temperature profiles and 
different heat transfer coefficients from those obtained bv 
assuming properties which are constant through the boundary 
layer. 

Assuming constant properties exist through the boundary 


layer, the boundary layer equations may be rewritten as: 





CONTINUITY 
Ju ۷ 
— + — = 0 (8) 
dX dy 
MOME NEUM 
130 تا‎ 10۳ 3 au 
ün 一 = - -一 :一 (一 -mv)) . (9) 
9 X dy pdx dy dy 
ENERGY 
ت3‎ ,+ + + + oa 
Ut ve = (at - v'T") ۳ 
9 X 3Y 0 3 ۱ ۱ 


The momentum and energy equations ((9) and (10)) are uncou- 


pled. This means that the momentum equation can be solved 
independently of the energy equation. In addition, the 
energy equation is linear. These two items simplify the 


process of solving the boundary layer equations. 

A number of existing methods to account for variable 
properties in boundary layers without film cooling will now 
be discussed. 


Mainstream gases consist of air and the combustion 


process bv-products (principally Carbon Dioxide, Carbon 
Monoxide, and water vapor). Viscosity and thermal conduc- 
tivity of this mixture can be determined though the semiem- 


pirical formula from WM EL 
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n nou: 
— Og 1 1 
E ee (11) 
i=l n i 
š m. 
j=1 J 1) 
where 
۱ )1 + Cre q T | 
ij BZ 
1) (8 + ۷ ۷ (۶ Eee 
1 J) 
This formula is for a mixture of n gases. Wi are the molec- 
ular weights and m; are the mole fractions. The exact same 


formula is recommended for thermal conductivity (u is 
replaced by k). Figure ll (White (1974)) illustrates the 


temperature variations of viscosity and thermal conductivity 


for these gases at low pressures. Notice that for a 200 R 
temperature variation, the variations in viscosity and 
thermal conductivity are not very dramatic. À power law of 
the form (Kays and Crawford (1980)): 
Hmix mix H didi 
—— € — 3 | حج‎ OUS 
U Do T 
mix, mix, e 
can be used to determine their temperature variations. For 
small temperature differences through boundary layers, say 


20 C or less, viscosity and thermal conductivity may each be 
considered to be constant. 

By considering the mainstream gases to be a perfect gas 
and from uniform static pressure through the boundary layer, 
the density variation becomes inversely proportional to the 


absolute temperature. 
D 工 
EAM LE (14) 


7 T 
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Comparsion of equations 13 and 14 shows that temperature 
variations result in larger density ratio changes than u or 


k changes. 
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III. PRELIMINARY CONSIDERATIONS 


The objectives in the thermal design of gas turbine 
blades are to obtain an acceptable stress level and design 
life. To do this, a detailed blade temperature distribution 
must be obtained. 

The thermal design of a cooled turbine blade is a 
complicated energy balance. The design process can be 
broken down into three major categories outlined below (in 
sequential order): 


(1) Predicting the heat flux from the hot mainstream 
gases to the blade s surface; 


(2) Performin a conduction 별 US for the detailed 
mapping of the blade s temperature distribution; 


(3) Predicting the heat flux from the blade to the 
internal cooling air. 


This thesis only addresses the first category: a relation 
is developed and verified between the two types of heat 


transfer coefficient used to predict wall temperature. 


۳ CONSTANT PROPERTY CONSIDERATIONS 
Equations 1 and 2 are used to calculate the blade's wall 
Cemperatume after the heat flux is obtained by solving the 


boundary layer equations. 


għ Be Tay Delos | (1) 


and 


qu ta DINE 2 T) c=) 
Both equations reflect the presence of film cooling. 

In developing thermal boundary layers with film cooling, 
an adiabatic wall boundary condition can be used to describe 
the convection process. Using this concept reduces the 


driving potential from a three temperature potential problem 
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to a two temperature potential problem. Equation l repre- 
Sents this concept. Here, he is an isoenergtic heat 
transfer coefficient. The adiabatic wall temperature, Toy 
is the wall temperature when there is no heat transfer to 
the wall. Adiabatic wall temperature is usually expressed 
nondimensionally in terms of the adiabatic film cooling 


effectiveness, n: 


T - T | 
aw ro 
Mo <a (150 
ې‎ 8 T 
re “ro 
Since the entire surface is considered adiabatic, ٢ is a 
global parameter. As a result, equation l is dependent on 


the entire convection process and the thermal boundary 
layer. 

Compared to the large temperature differences between 
the mainstream gases and the blade's surface, the surface 
temperature distribution can be considered relativelv isoth- 
ermal. An isothermal wall is ideal since thermal stresses 
are minimal. Equation 2 represents this concept. Here, h 
is defined by the freestream stagnation temperature and the 
wall temperature difference for a specific coolant tempera- 
ture. The coolant temperature is usually expressed nondi- 
mensionally as 8: 


Nc “To 
(16) 





0 ۲۷ 
W ro 


There has been considerable controversy over these two 
methods of predicting the wall's temperature. Two key 
issues of this controversy deal with how the temperature 


difference is defined and the effects of changing coolant 


parameters (temperature, density, mass flow rate, etc.). In 
equation و‎ when the temperature difference equals zero, 
there is no heat transfer. This is not necessarily true if 
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equation 2 is used. To account for viscous dissipation, the 
lunam stdagnatlon temperature, Ina. is replaced Dy the 
freestream recovery temperature, Tie The heat transfer 
coefficient, h, is then redefined as: 


٩1 ٤ e! (uio 


Changing the film cooling parameters will have different 
effects on these equations. Generally, the recovery temper- 
ature is determined from freestream conditions (temperature, 
velocity, density), however, if it is based on near wall 
conditions, the recovery temperature in a film cooled 
boundary layer is the same as the adiabatic wall tempera- 
Lone  Insequation ll; both he and Tj, are affected, whereas, 
sn equatieón 17 sly h is affected. 

For a given location, the heat flux will be the same for 
either method. As such, equation l is set equal to equation 
17 and a relation between the heat transfer coefficients 
ENS ES. Using the definitions of n and 8, and dividing by 
the heat transfer coefficient for the case of no film 


cocine, then yields: 


h ne 
— = —(1 - n9) (SN) 


This equation relates adiabatic and isothermal wall boundary 
conditions for properties that are constant through the film 
cooled boundary layer. 11715 2677 ج حط د5ہ‎ is plotted on 
eure 053۴. and Richards, 1978). Note the 
linearity between h/h, and 8. This linearity is due to the 
fact that the energv equation is linear for constant proper- 
ties. This equation is valid for a given downstream loca- 
tions where adiabatic conditions are represented by the 


horizontal intercept, and isothermal conditions by the line. 


Om Figure 12, the axis intercept points Nave SEN 
meaning. The h/h, axis Intercept Point 6 5615۹ nê 
transfer coefficient when the 09012717 تاد‎ o e eee 
8=0). The 8 


axis intercept point indicates a wall condition that is 


same temperature as the mainstream (Trc-Tr,, 


locally adiabatic. At this point, 8 equals the inverse of 
the effectiveness, ۰ 

Figure 12 is only good for one downstream location. AS 
you move further downstream, the slope of the curve becomes 
less and may eventually become horziontal indicating no 
cooling effect from the film. This is a result of a rising 
film temperature due to heat transfer, diffusion, and mixing 
with mainstream gases. As the film temperature rises, the 


value of © decreases and the value of N increases: 


B. -VARIABLE 807 1171 تظ‎ 05117151153 

With variable properties, the momentum and energy 
boundary layer equations become coupled because of the 
density term. in addikiom, the energy equation becomes 
highly G T Because of this, the linearity exhibited 
on Figure 12 no longer exists. Ligrani and Camci (1383) 
illustrated this nonlinearity that exists between h/hg and 
8. They also derived an empirical relation to relate the 
constant property data of Jabbari and Goldstein (1978) 988 
their variable property data. 

To find the relation between h/h, and 8 for the variable 
property case, we start with the defining equations for the 
heat transfer coefficients. 

For variable properties through the boundary layer, the 
variable property isoenergetic heat transfer coeffl- 
cient, há is defined using: 


Dm dv ço : 9 
٦ DU T (13) 


The heat fransier coefficient, h, is then defined using: 


S STO (20) 


The heat flux is now obtained by solving the variable prop- 
erty boundary layer equations. At a given location, the 
heat flux from equations 19 and 20 will be the same. 


Equating these gives: 


Dr 1 E 


٦ | 
aw Y co W 


f 
Now the adiabatic wall temperature is expressed in terms of 


an adiabatic film cooling effectiveness defined as: 


A 


A ㅠㅠ T 
ns حا‎ (22) 
R - T | 
Lc re 
٥ is defined as before (equation 16). By substitution and 


rearrangement the variable property relation between n/h, 


and 9 then becomes: 


一 =- Ñ(1 - n8) (23) 


h 
h he 
h 


O ho 


This equation is for one downstream location at one 


value of the density ratio, Dial? axis and one value of any of 


the following; blowing ratio,M; momentum flux ratio, I; or 
the velocity ratio, uj/u,. A different line results from 
equation 23 for each value of the effectiveness, n, and the 


normalized heat transfer coefficient, he/ho- The effective- 
ness, n, may have a strong dependence on the density ratio, 
کم‎ as indicated by the single row film cooling data of 
Pedersen, Eckect, and Goldstein (1977). 

Figure 13 is a domain map showing density ratio verses 
blowing ratio of the Pedersen, Eckect, and Goldstein (1977) 
data. Ihe box symbols indicate that effectiveness has 
little or no dependence on the density ratio for low blowing 
ratios. The circle svmbols indicate that effectiveness has 
a strong dependence on the density ratio. The triangles 
indicate that the effect of the density ratio on the effec- 


tiveness and heat transfer is highly nonlinear. 
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Figure 14 shows one possible effect of variable proper- 
ties on h/h, verses 9 data. Here, each of ehe 40066 linear 
curves represents equation 23 with a different n. The same 
hf /ho was used for each line to simplify Figure 14. Now, 
each line is also for one value of ۶۶06157۲۰ ti بت‎ Coe 
which is equivalent to T,/Tc. The variation of the density 


ratio with 8 is then given by: 





DEE 1 f. (36 

ACE = 
EE Te / To 9ت0‎ for ۰ 5 8 E The O fora given 
density ratio will then be a straight vertical line 
eminating from the horizontal axis. The intersection of 
this vertical line and equation 23 will be where a given 
h/h, data point is located. The curved line represents the 
collection SE the h/h, verses 8 data points. It curves 
downward because of the way the data collects. As the 


density ratio increases (vertical lines move to the right), 
the effectiveness increases.. This causes the O axis inter- 


cept point (1/n) to move to the left. 
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IV. RESULTS 


Using equation 23, the adiabatic effectiveness was used 
to deduce h/h, using the data of Pedersen, Eckect, and 
Goldstein (1977) for one row of injection holes. In addi- 


tion, the adiabatic effectiveness was deduced from h/h, and 
he/h, from the data of Ligrani and Camci (1985) for two rows 
of injection holes. 


In the following sections, the results are presented. 


B> DATA OF PEDERSEN, ECKECT, AND GOLDSTEIN (1977) 

Figure 15 shows the Pedersen, Eckect, and Goldstein 
11977) data used to deduce h/h,. Figure 13 gives the 
blowing ratio and the density ratio values for which n was 
selected for deducing h/h, - | 

A value of 1.0 was assumed for the h/h, axis intercept. 
Although in many cases, a reasonable assumption, Metzger, 
Carper, and Swank (1968) showed that this axis intercept may 
not be equal to 1.0. They showed that the injection has an 
effect on the heat transfer coefficient even when the 
coolant is injected at the freestream temperature, primarily 
due to boundary layer thicking. ۱ 

Normalized heat transfer coefficients as a function of 8 
are shown on Figures 16 thru 19 for the Pedersen,  Eckect, 
and Goldstein (1977) data. Four axes are used to denote 
four downstream locations (x/d). The light linear lines 
represent equation 23 for a different value of effective- 
ness, n. The dark lines represent the collection of h/h, 
verses O data points. 

ط٣ دا‎ TOT a blowing «ratio, M, of 0.213. This 
figure shows that /ط‎ verses 9 is steeper at small x/d and 
POEORESNNENAFEEC as. the distance from the injection site 
increases (increasing x/d). This figure shows reasonable 


linearity with some nonlinearity near the injection site. 
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Figure 17 is for a blowing ratio. Nous This 
figure shows the same trends as with Figure 16. Only mow, 
the nonlinearity near the injection site iS more pronounced 


and is also evident at a location further ۵9-15 dm 


٩157ه‎ ١د‎ 1 5 a DLONAME 67 ٥٠٣ OB AL" OD: Tois 
figure shows that Ene var 3 E h/h, with 8 is nonlinear. 
AS the distance from the injection site increases 
) 2: ۶2 5 2 ۳ ٥ h/h, variation with 8 becomes more 
hom onta Notice that, in general, at a particulas 
blowing Salon the slope becomes more negative as 8 
increases. 

Figure 19 is focla dbloninmetratio, TM, epe Jo vo This 


| figure shows a highly nonlinear behavior due to the 'lift- 
off" of film cooling jets from the wall, and possible pene- 


tration through the boundary layer. 


B. DATA OF LIGRANT SAND GSMGI (1985) 

Figure 20 (LiEramm and  Came+, 1985) illustrates the 
downward sloping trend for variable property nn verses ë8 
data. This figure also shows the linear h/h, verses 8 data 
(constant property case). Notice the similarity between the 
variable property trend here and che trends of the ۲ ۰ 7 ( 
Eckect, and Goldstein (1977) data on Figure 18. 

Using h/h_ and h,/h, from Ligrani and Camci (1985), n 
was deduced from equation 23. Results are plotted as a 


function ofewx/d in Figures 21 thru Bl, sand” as 2 1 رر‎ 





of E in Figuaes 82 theu 4% The later function, €, isa 
downstream distance normalized by the blowing ratio. L سيين‎ 
defined by: 
X e 1.909] U E ^ 
e A ~ (25) 
MS u 


© 


where S=nd/6 is the width of an equivalent two dimensional 
słot (Jabbari and Goldstein, 1978) . Goldstein (1971) gives 


two equations to -predict laterally averaged effective- 
ness, n. Por two dimensional falm Cooling. These 


predictions are given by: 


(26) سو جن ےر وت 
and‏ 

= 1.9pr°/> 

"1 + 0.329g0-38 ا‎ 
where | 

pt 1.5(1077) (Re. (-E))sinr (23) 


U 


So 


These predictions are plotted on the n verses £ graphs 
(Figures 32 thru 42). 

Figures 21 thru 42 are for different blowing ratios (M = 
17242, 1.003, 0.481, 0.661, 1.557), momentum flux ratios (I 
5-۹-2970 wand velocity ratios (u_/u, = 0.630, 


IRS 7722.). In general, these figures show that as the 
blowing ratio, momentum flux ratio, and velocity ratio 
increase, the adiabatic film cooling effectiveness 
increases. The largest increase in n occurs for the highest 
value of O. This corresponds to the case of the highest 
density ration. On. Figures 21 thru 31, the crossover point 


between the Ol and 02 curves moves downstream (from x/d = 20 


to x/d = 40) as the momentum flux ratio and velocity ratio 


increase. On Figures 32 thru 42, as the blowing ratio, 
momentum flux ratio, and the velocity ratio decrease, the 
data moves to a higher value and spreads out. All data 


points on the n verses & graphs are below the predictions of 
equations 26 and 27 but do follow the same trend as the 


predictions. 
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V. SUMMAR 


A relation was derived that relates the normalized heat 


transfer coefficients h and he. These heat transfer coeffi- 
cients are for isothermal wall boundary conditions, h, and 
adiabatic wall boundary conditions. hg. Each are for when 


variable properties are considered to exist through the 


boundary layer. This relation was presented as: 
h he ^. 

— = —(1 - ne) (23) 
ho ho 


When the effects of variable properties are important in 


film cooled turbulent boundary layers, the variable property 


effectiveness, n, and the above equation are valid for one 
density ratio and blowing ratio only. 

The normalized heat transfer coefficient, h/ho, Was 
deduced from the Pedersen, Eckect, and Goldstein (1977) data 
on injection with one row of holes. The adiabatic film 
cooling effectiveness, n, was deduced from the Ligrani and 
Camci (1985) data on injection with two rows of holes. The 


data deduced from equation 23 showed the nonlinear trend of 


Varied. 


oo > 


h/h, verses O0 and as the density ratio, oc /P 

The effect of density ratio on h/hg verses © results in 
a change of the characteristic curve from linear to nonli- 
near. This change is not always evident. It's existance or 
extent is dependent upon the injection parameters (blowing 
ratio, momentum flux ratio, and velocity ratio) and the 
downstream distance. This effect is due to the nonlinearity 
introduced into the energy equation and the coupliing of the 
momentum and energy equations when variable properties 


through the boundary layer are considered. 
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VI. RECOMMENDATIONS 


Experiments need to be conducted to determine the accu- 
racy and dependabilitv of the proposed variable property 


relation on injection system geometry, injection flow rate 


parameters, fluid mechanics of the flow, the thermal 
boundary layer conditions, and other parameters. In order 
to completely verify the proposed equation, Simultaneous 


measurements of variable property effectiveness and variable 
property heat transfer coefficients (based on the wall and 
recovery freestream temperature difference) need to be 


conducted. 
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